Fluctuation theorems: Work is not an observable 
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The characteristic function of the work performed by an external time-dependent force on a Hamil- 
tonian quantum system is identified with the time-ordered correlation function of the exponentiated 
system's Hamiltonian. A similar expression is obtained for the averaged exponential work which is 
related to the free energy difference of equilibrium systems by the Jarzynski work theorem. 
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Deep relations between nonequilibrium fluctuations 
and thermal equilibrium properties of small systems have 
recently been discovered and formulated in terms of so- 
called fluctuation theorems [1, 2]. These theorems are 
not only of basic theoretical relevance but provide novel 
ground for experimental investigations of small systems 
in physics, chemistry and biology [3] 

In this rapid communication we want to preclude a pos- 
sible confusion about the notion of work in the context 
of fluctuation theorems for quantum mechanical Hamil- 
tonian systems. To be precise, we consider a quantum 
system which was in thermal contact with a heat bath at 
inverse temperature (3 until a time to- Then the contact 
with the bath is switched off and a classical force acts on 
the otherwise isolated Hamiltonian system until the time 
tf according to a prescribed protocol. We demonstrate 
that the exponential average of the total work performed 
on the system as well as the characteristic function of 
this work are given by time- ordered correlation functions 
of the exponentiated Hamiltonian rather than by expec- 
tation values of an operator representing the work as a 
pretended observable. 

For a system that evolves under the exclusive influence 
of a time dependent Hamiltonian H(t) from an initial 
thermal equilibrium state 



p(0) = Z(0)- 1 exp{-/3 J ff(0)}, 



(1) 



at time to = until a final time t = tf, the work per- 
formed on the system is a randomly distributed quantity 
w. Its statistical properties follow from a probability den- 
sity p(w) or, equivalently, from the corresponding char- 
acteristic function G(u) which is defined as the Fourier 
transform of the probability density, i.e. 

G(u) = j dwe mw p(w) (2) 

We will demonstrate that the characteristic function is 
given by the following quantum correlation function: 

G(u) = ( e »ff(*/) e -»ff(o)} 

ee Tie mHH( - t f ) e-' luH ^p(0) 

(3) 

= Tr e iuH MU(t f )e- iuH ^ p(0)U + (t f ) 

= Tie iuH ^l7(t/)e-( <u+ ^ K Wl7 + (t / )/Z(0) 



where Tr denotes the trace over the system's Hilbert 
space TL, U(t) the unitary time evolution governed by 
the Schrodinger equation ihdU(t)/dt — H(t)U(t) with 
U(0) = 1, and H H (t) = U+{t)H{t)U{t) is the Hamilto- 
nian in the Heisenberg picture. The third equality follows 
from the second line by the cyclic invariance of the trace 
and the last line follows with eq. (1). 

The characteristic function has the form of a time- 
ordered correlation function of the two operators 
exp{iuH(tf)} and exp{— iuH(0)}. We note that this 
correlation function in general differs from the aver- 
aged exponential of the difference of the Hamiltonians 
W = Hji(tf) — H(0), which sometimes is referred to 
as the operator of work [4]. It is possible though to 
formally rewrite the characteristic function in terms of 
the difference between the Hamiltonians. In the second 
line of eq.(3) the product of the operators exp{iuiJ#(i)} 
and exp{— iuH(0)} occurs in chronological order and 
may be written as exp{iuHn (t)} exp{— iuH(0)} = 
T> exp{iuH n{t)} exp{— iuH(0)}. Under the protection 
of the time ordering operator T > the usual rule for expo- 
nentials of commutative quantities holds [5] to yield the 
following equivalent forms of the characteristic function 
of work 



G{u) = Trr >e m ^(*^- ff ( o »p(0) 
C tf dH H (s) 



Tr T> exp{iit / 



ds 



ds}p(0) (4) 



The second equality is a consequence of the known fact 
that the total derivative of the Hamiltonian in the Heisen- 
berg picture coincides with its partial time derivative. 

The averaged exponentiated work (exp{— [3w}} is ob- 
tained from the characteristic function by putting u — i[3, 
cf. eq. (2). Using the correlation function expression (3) 
together with the canonical initial density matrix (1) we 
immediately recover the Jarzynski equation in its known 
form [6] 



Z(tf) 
Z(Q) 



(5) 



where Z(tf) — Tr e~P H ( tf ' is the partition function of a 
hypothetical system with Hamiltonian H(tf) in a Gibbs 
state at inverse temperature (3. 



By replacing the quantum correlation function by the 
corresponding correlation function of a classical Hamilto- 
nian system the characteristic function of the work per- 
formed on the classical system is obtained. Its inverse 
Fourier transform yields the known classical expression 
for the probability density of work, cf. Ref. [7] 

PdW - Z^iO) J dz(0) 

5{w-[H(z(t f ),t f )-H(z(0),0)\) 

X(= -/3ff(z(0),0) (6) 

where Z C \(Q) = J dzcxp{— /?i/(z, 0)} denotes the classi- 
cal partition function, z = (p, q) a point in phase space 
which serves as the initial condition of the trajectory z(i) 
evolving according to Hamilton's equations of motion. 

The fluctuation theorem has long been known for a 
sudden switch of the Hamiltonian of a classical system 
[8]. For a quantum system with a Hamiltonian chang- 
ing from H at time i = 0~ to Hi at tj = + the 
time-evolution operator becomes U(tf) = 1. The char- 
acteristic function (3) then simplifies to 

G(u) = Tr e tuIil e~ iuHo e^ /Z(0) (7) 

and the Jarzynski equation (5) again follows with u = i(3. 
In all nontrivial cases, when the two Hamiltonians do 
not commute, the averaged exponential of the difference 
operator Hi — Ho docs not yield this result. 

The proof of eq. (3) essentially follows an argument 
given by Kurchan [9], see also Refs. [10-12]. It is based 
on the elementary observation that two energy measure- 
ments are required in order to determine the work per- 
formed on the system by an external force. In the first 
measurement, the energy is determined in the initial 
Gibbs state Z(O)" 1 exp{-/3i/(0)}. The outcome of this 
measurement is one of the eigenvalues e„(0) of the Hamil- 
tonian H(0) with the probability 

p n = exp{-/3e„(0)}/Z(0). (8) 

After the measurement the system is found in the 
corresponding eigenstate c/?„(0) of H(0) satisfying 
H(0)ip n (0) = e„(0)</?„(0). This state evolves according to 
il>(t) = U(t)ip n (0) until the second energy measurement 
is performed at the time tf. It produces an eigenvalue 
e m(tf) with the probability 

p(m,t f \n) = \(<p m (t f )\U(t f )<p n (0))\ 2 , (9) 

where (•[•) denotes the scalar product of the Hilbert space 
H. Here, e m (tf) and ip m (tf) are the eigenvalues and 
eigenfunctions, respectively, of the Hamiltonian H(tf). 
Hence, the energies e m (tf) and e n (0) are measured with 
the probability p{m, t\n)p n such that the probability den- 
sity of the work, which is the difference of the measured 
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energies, becomes 
p( w ) = ^2 S ( W - [ e m(tf) - e„(0)]) p(m,t f \n)p n (10) 

One then finds for the characteristic function from the 
definition (2) 

G(u) = ^e m ( e "^- e "(°» 

n.m 

x (<p m (t f )\U(tf)<p n (0)) 
x (M0)\U + {t f )<p m (t f ))e-^/Z(0) 
= J2^m(t f )\U(t f )e-^p(0) Vn (0)) 

n.7n 

x (^ n (0)\U + (t f )e iuH ^^ m (t f )) 
=TiU{t J )e- iuH{ - 0) p{Q)U + {t J )e iuH{ - t ' ) 

In going to the last line, the sum over the complete set 
of eigenstates {ip m (t)} was written as the trace and the 
completeness of the eigenstates <£„(0) was used. By use of 
the cyclic invariance of the trace the quantum correlation 
function expression (3) for the characteristic function is 
proved. 

This expression for the characteristic function contains 
all available statistical information about the work per- 
formed by an external force on an isolated quantum sys- 
tem, such as the averaged exponentiated work, cf. eq. 
(5). All moments of the work follow from the derivatives 
of the characteristic function. 
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